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Abstract 

CN ■ We study the effect of parameter uncertainty on a stochastic diffusion modei, in par- 

ticular the impact on the pricing of contingent claims, using methods from the theory of 
Dirichlet forms. We apply these techniques to hedging procedures in order to compute the 
sensitivity of SDE trajectories with respect to parameter perturbations. We show that this 
jrt ■ analysis can justify endogenously the presence of a bid-ask spread on the option prices. We 

also prove that if the stochastic differential equation admits a closed form representation 
^^ ■ then the sensitivities have closed form representations. 

We examine the case of log-normal diffusion and we show that this framework leads to 
a smiled implied volatility surface coherent with historical data. 
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Q_| , Theory using Dirichlet Forms, Bid- Ask Spread, Value- at- Risk. 
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1 Introduction 

^H \ 

The purpose of this article is to study the effect of parameter uncertainty on asset pricing. 

\q • If we consider, for instance, the Black and Scholes model, the price of a stock is modeled as 

^sO ! a log-normal process. This diffusion process depends on a random source, represented by a 

Brownian motion, that models the market uncertainty and the unknown future evolution of 
the stock price. This diffusion depends also on two parameters, usually known as the drift 
and the volatility. The latter plays a crucial role in asset pricing. As such, a second source 

r—i \ of uncertainty appears since the values of these parameters have to be estimated by financial 

J> ■ agents. In practice, in order to sell an option, one has to fix a price using his pricing method, 

which depends on the diffusion parameters. Since these diffusion parameter are unknown, one 
has to plug into his pricing model their estimated values. However, these estimated values are 
uncertain, in other words, they are random variables. A natural question arising is therefore 
the impact of this second source of uncertainty. 

A large literature exists on this subject. We may in particular refer to the uncertain volatility 
model (UVM) which take into account the difficulties of calibrating the volatility in the Black 
and Scholes model, see Avellaneda et al. [2] and Lyons [25] • In UVM, the parameter, i.e. the 
volatility, is an unknown parameter which belongs to a interval. However, this kind of approach 
neglects the way practitioners implement their pricing method by using estimated parameters. 
On the other hand, practitioners often overlook the inherent uncertainty in their estimated 
parameters. 

In our model, we assume that the stock price follows a stochastic differential equation (SDE) 
with fixed but unknown parameters. As a consequence, only one random source exists in the 
market, which is therefore complete. However, the financial agents, assumed to be option sellers 
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in our analysis, know that the right model is an SDE but they do not known the true values 
of the parameters. Therefore, they estimate the parameters and use these estimates to price 
the options. In other words, the parameters are replaced by statistical estimators of the true 
parameter values, that is, they are random variables characterized by a variance and a bias. As 
a consequence, the option seller uses a different SDE, i.e. an SDE with different parameters to 
compute the option prices. It follows that the probability space used in this paper must describe 
two random sources: the market uncertainty, i.e. the Brownian motion driving the SDE, and 
the parameter uncertainty, that affects the calculations made by the option seller. However, 
the two uncertainties do not play the same role. Since the law of the market uncertainty is well 
defined while the law of the parameter estimators is poorly known. Because of this, we assume 
known only the biases and the variances of these estimators. Our goal is to study how the bias 
and the variance of the parameter estimators affect the basics of the asset prices and the related 
option prices. 

An important difference between our work and the UVM model is that we consider a large 
class of underlying diffusions and not uniquely the log-normal one. We assume that the true 
underlying diffusion is unique and follows a stochastic differential equation with continuous 
paths, 

dS t = S t lidt + S t <j(t, S t ) dB t . 

This diffusion can depend on several parameters, and not only on the Black-Scholes volatility, 
that are to be estimated by the option seller. We assume that the uncertainty is on the volatility 
and carried by the parameters used to describe this volatility. Our method leads to a option 
pricing calculus that enables the option seller to control the risk due to parameter uncertainty 
assuming an acceptable residual risk that can be evaluated using a Value-at-Risk technique, we 
propose a principle to price an option and we find a selling price such that the probability of 
loss of the option seller is smaller that a given tolerance a, we explain our pricing principle in 
Section 3.2. 

In the UVM case, the strategy is to apply stochastic control techniques to super-hedge a 
contingent claim. The principal drawback of the super-replication technique is that the super- 
hedging cost is too high and the corresponding strategy too conservative. As an example, 
Kramkov [23] shows that the super-hedging cost of a call corresponds to the value of the option 
under the least favorable martingale measure. Likewise, Bellamy and Jeanblanc [3] prove that 
in a jump-diffusion model the price range for a call option corresponds to the interval given 
by the no-arbitrage conditions that is the super-hedging strategy of a call option is to buy 
the underlying. It seems therefore that super-hedging is not an effective methodology since 
the option buyer will prefer to buy the underlying rather than pay the same price to have an 
option. Many attempts to overcome this problem have been studied in literature. Avellaneda et 
al. [2] constrains the volatility between two levels. Another approach is the pricing via utility 
maximization, see for instance El Karoui and Rouge |14| . or via local risk minimization, see 
for instance Follmer and Schweizer |15j . We also recall the paper of Denis and Martini |llj . 
where a more general framework is introduced in order to study super-replication and model 
uncertainty. Their methodology makes use of Choquet capacities and is related to our approach. 
Recently, a generalization of this approach is proposed by Soner et al |27j . 

In our analysis, we assume that the option seller has to take some risk in order to propose 
a competitive price. In other words, he has to take the risk of losing money, with a small 
probability, in order to have the opportunity to make money with the contract (selling the 
option and hedging it). We quantify this risk thanks to a Value-at-Risk technique. In a similar 
way as in UVM, we find distinct buying and selling prices due to the asymmetry of risks between 
buyer and seller. 

In practice, the risk assessment is often carried out in terms of sensitivity with respect to a 
deterministic variation of model parameters. These sensitivities are usually called greeks and op- 
tion sellers try to manage their portfolios in order to control the risk related to these sensitivities 
using risk measures like historical Value-at-Risk. However, this approach leads to mathematical 



difficulties, for instance in dealing with infinite dimensions, and practical difficulties. 

An alternative method that we will use here is based on potential theory (see Albeverio 
[I], Bouleau and Hirsch |5] and Fukushima et al. |17j ) suggested by Bouleau [6]. We assume 
that the uncertainties are small random variables added to the true values of the parameters. 
The laws of these estimators are generally poorly known, so we assume that only the two first 
moments are known, i.e. the biases and the variances/covariances. When the uncertainty on 
a parameter is small, it is reasonable to neglect high moments. This is also justified by the 
fact that parameter estimation is often difficult owing to the non-linearity, the complexity of 
the equations and to the poor quality or the shortage of data. The theory of Dirichlet forms 
leads to an "error theory" allowing the computation of the two first moments of the random 
variable that describes the option seller's profit or lost due to the parameter uncertainty. This 
theory combines the variance to the "carre du champ" operator and the bias to the generator 
associated to the Dirichlet form, see Bouleau [7]. 

To summarize, the main result of this paper is to give closed forms for the corrections on 
the price of contingent claims due to the uncertainties on the parameters that are modeled 
as statistical estimators, i.e. random variables, with small variances and biases. These are 
corrections for an investor seeking to compensate the bias on the option price and accepting 
a residual risk, both due to the uncertainty on parameters. We define the risk taken by the 
option seller in a statistical way in accord with a VaR risk measure. A direct consequence of 
this pricing method is the separation of the buy and sell prices due to the asymmetry between 
the buyer and the seller with respect to the risk on parameters. A surprising result is that a 
systematic bias exists with respect to the pricing without uncertainty even if the parameters 
are unbiased. Indeed, this bias is due to the non-linearity of the payoff. Our analysis covers a 
large class of stochastic process with continuous paths and, in particular, equity models with 
local or stochastic volatility. 

Finally, we present an example based on the Black-Scholes model, with uncertainty on the 
volatility parameter. We show that that our analysis leads to a smiled implied volatility surface 
and we give the explicit formula for vanilla option prices and the bid-ask spread. 

The paper is organized as follows. In Section 2, we resume the notation used in the paper. 
In Section 3, we briefly recall the error theory based on Dirichlet forms techniques. Section 
4 deals with the problem of uncertainty estimation while Section 5 is devoted to the study of 
the impact of uncertainty on the diffusion model. We analyze the profit and loss process and 
compute its law depending on both the underlying diffusion and the parameter uncertainty. In 
Section 6, we introduce a pricing principle to over- hedge the contingent claim and exhibit the 
bid and ask prices. In Section 7, we give an example based on log-normal diffusion without 
drift. We exhibit the bid and the ask prices and we prove, under certain hypotheses, that 
the implied volatility exhibits a smile behavior. Section 8 presents a second detailed example 
based on constant elasticity of variance model in order to show the powerful of the method in 
a practical case. 

2 Notation 

This section fix the notation used throughout the paper. We consider two probability spaces 
(fil, J- 1 , {J't}t£[o,T]^i) an d (fi,.F, P). The reference filtration {J 7 / } ie [o^i of the first probability 
space (fiij.F 1 , {JF^ }te[0,T]>Pl) satisfies the usual condition and is generated by a Brownian 
motion B. In this space, we define a continuous time financial market model consisting of a 
risk- free asset, whose price process is assumed for simplicity to be equal to 1 at any time, and a 
risky asset (or stock) of price process S = (St)t>o- We will denote by Ei the expectation under 

The second probability space (£l,J~, P) is the sample space used to define the parameter 
estimators. We will denote by small letters, e.g a, the true parameters of the model; by capital 
letters, e.g. A, the estimators of these parameters and by capped capital letters, e.g. A, the 



estimated values of these parameters. But, to avoid any confusion, the true volatility value will 
be denoted by <r, its estimator by ? and its estimated value by <f. As a consequence, A will be 
a random variable defined in the space (Q,J-,¥), when A will be a fixed value resulting from a 
particular observed dataset. 

We consider the product space ($7, J 7 , P) = (Oi xfl,/ 1 ® .F,Pi X P). We assume that the 
Brownian motion B remains a Brownian motion under the probability space (£1, F, P) and that 
all estimators defined in the probability space (Q, F, P) are independent with respect to the 
filtration generated by the Brownian motion. 

We need to introduce an auxiliary probability space denoted by (Q, J-, P). This space is 
a copy of estimator probability space (O, F, P). We will denote by E the expectation under 
the probability P. This auxiliary probability space is introduced only to perform mathematical 
proofs and has no influence on final results. 

Finally, we introduce a representation basis, denoted by {4>i(t, x)}i£j\f, of the space of volatil- 
ity function. All functions (j>i belong to C 1 ' 2 , {4>i}t£^ and their derivatives are Lipschitz by a 
same constant and square integrable in x. The space of admissible volatility functions, denoted 
by L 2 a , is 

L 2 = < f(t,x) = Y^bi<f>i(t,x) such that (& 1? . . . ,b n , . . .) G £ 2 and f(t,x) > £ > OV(i,x) > 

We equip the space with the distance d 2 induced by the sequence representation. 

3 Mathematical tools 

We begin by giving a general introduction to the study of sensitivity with respect to a stochastic 
perturbation and a formal definition of the framework that we will use, i.e. error theory based 
on Dirichlet forms, in accord with Bouleau [6]. In this survey, we follow [7]. 

We consider a function F(ux, U2, •••, u n ) depending on the parameters u = (u±, U2, ■■■, u n ) 
that we assume approximately known with some uncertainty. That is, we assume that we have 
a statistical estimator U = (U\, U2, — , U n ) of the vector u. We consider that the law of U is 
unknown and we suppose known only the variance and bias. We assume that the function F is 
sufficiently regular and we seek to evaluate the impact on F of the uncertainties on u. The first 
study of this problem goes back to Gauss who proved the following expansion for the variance 
of F(U) if the uncertainties are small compared with the parameters values: 

Vax[F(U)] =J2^(u) J^(«) Covar[^, U d ] 

However, this relation is valid only if the number of parameters is fixed and if F has an explicit 
formula, in particular, it cannot be defined via a limit. Another main problem is that we do 
not know u but only U, i.e. the estimator of u. As a consequence, we cannot evaluate the 
derivatives of F at u. Some natural questions that arise are: can we use the same approach if 
we know only the random variable Ul That is, can we prove that 

Var[F(C/)] = £ H(EO |£V) Covar[f/ 4 , U s \ ? 
i,3 J 

Is it this relation true when F is defined implicitly, for instance by a L 2 -limit? Can we study the 
bias of F(U) caused by a non- linear function Fl The classical Gauss theory does not answer to 
all of these problems. Therefore, we turn us to the error theory based on Dirichlet forms in order 
to go beyond these problems. In particular, we need to evaluate the effect of a stochastic error 
on more complex objects such as stochastic integrals. For this, we consider an estimated value 



U of a parameter u with a small uncertainty Y, on which we compute a non-linear function F. 
Our estimate U is replaced by the estimator U = U + Y, i.e. a random variable with a mean 
squared error MSE[Y] and bias Bias[Y]. Appling Taylor's expansion to a function F, we find 

MSE[F(J7)] = I \{F{U) - F(U)} 2 ] = \f'(U)\ MSE[U] + higher orders 
Biaa[F(U)] = E[F(U) - F(U)} = F'(U) Bias[U] + -F"(U) MSE[U] + higher orders . 

Clearly, the mean squared error MSE[i ? (C7)] dominates the variance V&r[F(U)] defined as 

E {F(U) — K[F(U)]} 2 . For the centered random variable U the mean squared error and 

the variance take the same value. Moreover, the difference MSE[.F([/)] — V&t[F(U)] is equal to 
the square of the bias, that is, it is negligible in our expansion. As a consequence, we will use 
the mean squared error as a conservative estimation of the variance. 

If we suppose the uncertainty small, we can truncate the Taylor's expansion and find two 
closed transport formulas for the bias and the mean squared error. These are known in literature, 
since the mean squared error has the same transport formula as a "carre du champ" operator on 
a probability space equipped with a local Dirichlet form, while the bias has the same transport 
formula as the generator of the semigroup associated to the Dirichlet form, see for instance 
Bouleau and Hirsch [5|. The main advantage of this comparison is that the carre du champ 
operator and the generator of the semigroup are closed operators with respect to the graph 
norm: see for instance Fukushima et al. |17| . Therefore, they are good operators for studying 
objects defined by limits, such as stochastic integrals. 

The axiomatization of this idea was introduced by Bouleau [6] as follows: He defined an 
error structure as a probability space equipped with a local Dirichlet form with a carre du 
champ operator. 

Definition 3.1 (Error structure) 

An error structure is a term ( Q,, F, P, ID, F I , where 



Q,, J-, PI is a probability space; 

• B is a dense sub-vector space of L 2 (Q,, J 7 , PI; 

• r is a positive symmetric bilinear function from D x D into L 1 I £7, J 7 , PI satisfying the 

functional calculus of class C 1 ClLip, i.e. if F and G are of class C l and globally Lipschitz, 
and U and V belong to B, then F(U) and G(V) belong to B and 

(1) T [F(U), G(V)} = F'(U) G'{V) T[U, V] P a.s.; 

• the bilinear form £[U, V] = iE \T[U, V]] is closed; 

• the constant 1 belongs to B and £[1,1] =0 

We generally write T[U] for T[U, U] and we denote by T [[/]([/) the particular realization 
of the random variable T[U] corresponding to observed data, i.e. r[C/](C/) is the conditional 
variance of the estimator U given the particular dataset used to estimate it and such that the 
estimated value of u is equal to U. With this definition, £ is a Dirichlet form and V is the 
associated carre du champ operator. The Hille-Yosida theorem guarantees that there exists 
a semigroup and a generator A that are coherent with the Dirichlet form £, see for instance 
Albeverio [1] and Fukushima et al. |17| . This generator A : T>A — > L X (P) is a self-adjoint 



operator, its domain T>A is included into D, and this operator satisfies, for F € C 2 , [/ € 2?^4 

and V[U] G L 2 (P): 

(2) ^ [F([/)] = F'(C7) ^[C/] + ^F"{U) T[U] P a.*. . 

Moreover, it is a closed operator with respect to the graph norm. We denote by «4[{7]({7) the 
particular realization of the random variable A[U] corresponding to observed data, i.e. ^4[C/](C/) 
is the bias of the estimator U given the particular dataset used to estimate it and such that the 
estimated value of u is equal to U. We emphasize two important points related to this theory. 
First of all, error structures have nice properties: in particular, it is possible to prove that 
the product of two or countably many error structures is an error structure (see Bouleau [7]). 
Furthermore, the elements of the space B are square integrable random variables. Therefore, this 
theory is coherent with the stochastic nature of the estimators U, and we consider the operator 
A and r as the conditional bias and mean squared error given the value of the estimator U: for 
details we refer to Bouleau [6], chapters III and V, and [7J. 

The classical example of an error structure is the Ornstein-Uhlenbeck one: 

Definition 3.2 (Ornstein-Uhlenbeck error structure) 

The Ornstein-Uhlenbeck structure is (R, £>(R), G, H , U — > (U ! ) I, where G represent a 

unidimensional Gaussian law, H 1 is the Sobolev space of order 1 and the carre du champ operator 
is just the square of the first derivative. 

The main drawback of the carre du champ operator is its bi-linearity, which makes com- 
putations awkward to perform. An easy way to overcome this drawback is to introduce a new 
operator, called the sharp operator. This operator is a specified choice for the gradient in the 
terminology of Bouleau and Hirsch (see Bouleau and Hirsch |3], section II. 6.). We recall the 
definition of the sharp operator associated with V. 

Proposition 3.1 (Sharp operator) 

Let (Cl, JF, P, B, r) be an error structure and let (O, J 7 , P) be a copy of the probability 
space (Q, J 7 , PI. We assume that the space B is separable. Then there exists a linear operator, 
called sharp and denoted by ( )# : B — > L 2 (P x P) ; with the following two properties: 

• \/U € B, r[£7] = E (U&) , where E denotes the expectation under the probability P; 

. MU € W and F £ C 1 n Lip, (F(U U ... , U n )f = YTi=i (jg ° u) U* . 

The sharp operator is a useful tool when computing T because it is linear, whereas the carre 
du champ operator is bilinear. Moreover, the sharp operator is closed, then, for instance, we 
can exchange the sharp operator and the integral sign. The proof of this fact proceeds by an 
approximation of the integral by a sum, then we apply the sharp operator and finally we take 
the limit using the closeness of the sharp operator, see for instance Bouleau jU] section VI. 2. 

We give a detailed example of the computations of sharp, carre du champ and bias operator 
in Section 8. To make our approach more understandable for a non specialist of Dirichlet forms, 
we explain in this example the computation of each operator. The classical approach is to 
compute the sharp operator. This operator is just an auxiliary step needed to the computation 
of the carre du champ, that represents the variance in error theory framework. Finally, the 
knowledge of variance operator T allows us to compute the action of the bias operator A. 

We observe that the error theory based on Dirichlet forms restricts its analysis to the study 
of the first two orders of error propagation, i.e. the bias and the variance. This fact is justi- 
fied by the lack of information on the parameter uncertainties, generally given by the Fischer 



information matrix, that is often quite limited. The study of higher orders is a very difficult 
problem for both mathematical and practical reasons. From the mathematical point of view, 
it would be necessary to study chain rules of higher orders, involving skewness and kurtosis, 
and to prove that the related operators are closed in a suitable space. However, the crucial 
problem remains to have sufficiently accurate estimates for the higher order uncertainties. This 
statistical obstacle cannot be overcome easily. Therefore, we decide to restrict our study to the 
two first orders. 

We conclude this survey recalling two direct consequences of this approach. 

Corollary 3.2 (Impact of uncertainty) 

If we neglect moments higher than the second one, i. e. we approximate all random variables 
by the Gaussian one with same expectation and variance, the impact of uncertainty on the 
parameter u transforms F{U) into a Gaussian random variable of the form 

(3) F(U) A F(U)+A[F(U)](U) + y/r[F(U)](U) G, 

where G is a standard Gaussian variable. 

If we want to be conservative, in particular if we suppose that the Gaussian approximation 
is not valid, then we can apply left-hand tailed Chebyshev's inequality (also known as Cantelli 
inequality), which leads to the following statement. 

Corollary 3.3 (left-hand tailed Chebyshev's inequality) 

The random variable F(U) verifies the following inequality for all k > 1. 

1 



(4) 



F(U) - F(U) - A[F{U)]{U) > k Jr[F(U)](U) 



< 



1 + k 2 



4 Estimation, calibration and uncertainty 

In this section, we discuss briefly the classical methods to determine the model parameters and 
how to evaluate their uncertainties. 

Statistical estimation. The classical approach deals with estimating the parameters using 
historical data of the stock price. The goal of this method is to construct an estimator, that 
takes the historical data as input and gives as result an estimated value for each parameter. 
Then, an estimator is a function mapping the historical data probability space (sample space 
throughout the sequel) into the set of possible estimated outcomes. It is important to remark 
that an estimator is a random variable whereas an estimated value is a constant. We can 
associate different attributes with an estimator, like its bias, its variance and its mean squared 
error. We recall that all these attributes depend both on the estimator, i.e. the rule to find 
the estimated value, and on the sample, i.e. the historical data. In this sense, maximum 
likelihood estimation plays a central role between statistical estimation methods. We recall 
the asymptotically efficiency between their properties, meaning that the maximum likelihood 
estimator reaches the Cramer-Rao bound, that is a lower bound on the variance of a parameter 
estimator. Bouleau and Chorro [8] have proved that the operator T is the inverse of the Fisher 
information matrix, i.e. the inverse of Cramer-Rao bound, in the finite dimensional case. 

Calibration on liquid options. The main drawback of statistical estimation is that the diffu- 
sion process depending on the estimated parameters is generally not consistent with the market 
prices of vanilla options. But the liquidity of these options is large enough to consider that their 
prices represent equilibrium prices. As a consequence, a financial model has to be consistent 
with these prices. The method of calibration is to find which sets of parameters are consistent 
with market prices of vanilla options. A general result due to Dupire |13| shows for instance that 



a unique local volatility model is coherent with the surface of option market prices across stikes 
and maturities. However, Dupire formula is an ill posed problem. A regularization method is 
then required in order to stabilize the solution. The classical regularization method is to add 
a penalization term and to find a tradeoff between consistency with data and smoothness, see 
for instance Cont and Tankov |10] , 

Calibration method gives a rule to find the best set of parameters, i.e. the set of parameters 
that minimizes a given optimization principle. This rule is a function from the sample space 
into the possible outcomes, like the estimator in the case of statistical estimation. But it is 
difficult to adapt the previous methodology to calibration case, since the calibration function 
from the sample space into the possible outcomes is not explicit. 

A different approach can be however proposed. It is well-known that, even if option volumes 
are grown exponentially in the last thirteen years, their relative spread remains large compared 
with assets relative spread. The relative spread is defined as the ratio between the bid-ask 
spread and the mid price. As a consequence the equilibrium price of a vanilla option is known 
with an uncertainty, measured by the bid-ask spread. This uncertainty has to play a role in the 
calibration methodology. 

We propose then an easy method to transfer these price uncertainties into an uncertainty on 
the calibrated parameters. The first step is to perform the calibration methodology neglecting 
the bid-ask spread, i.e. for any option we fix the price at the mid-price. Afterward, we fix an 
option j on the basket used to calibrate, we shift its price to ask price and we recalibrate. The 
new set of calibrated parameters represents a stress of the previous set. We compute then the 
difference between the two sets and we interpret it as a standard deviation of the calibrated 
parameters submitted to a stress into the price of option j. So we can reconstruct the variance- 
covariance matrix with respect to the random source in option j. We perform this procedure 
for all options into the basket and we can reconstruct the global variance-covariance matrix. 

This method can be easily performed given that only a rough approximation of the stressed 
calibration set is needed. Moreover, the initial calibration set can be used as a starting set 
to perform stressed calibrations. A corollary result of this approach is to evaluate the fit of a 
model to the real data, since a large variance-covariance is the sign of an overfitting situation 
due to a too large number of parameters or to an unconformability of the model behavior to 
real data shape. 

Pure calibration method. Finally, a third method to define parameters uncertainty in our 
framework is via a pure calibration approach. In other word, we can study theoretically the 
behavior of the model depending both on the estimated parameters and their uncertainty. 
Given the pricing formula or the pricing rule, for instance the simulation using Monte Carlo, it 
is possible to study the parameters sets coherent with vanilla options. 

The crucial remark of this section is that statistical estimation and calibration produce a set 
of estimated parameters, that is the best approximation of true parameters given all historical 
data. As a matter of fact, these best estimated values are both 

Fixed Values since historical data are fixed and known. Then the rule to estimate, via sta- 
tistical estimation or calibration, gives us a numerical result. 

Random Variables since historical data are in fact a given realization. Then the dependence 
with respect to the simple space cannot be neglected. 

As a consequence, in the pricing methodology the estimated parameters has to be considered 
fixed and this set of parameters is the only one that can be used to perform all computation. 
However, an estimation risk appears and the dependency on the given realization of historical 
data becomes the kernel question to be analyzed. The estimated volatility is then a fixed value 
(denoted by q) on the probability space (jli, J 71 , Pi) and a random variable (denoted by <j) on 
the probability space (O, JF, P). 



The main objective of the present paper is to investigate the impact of this double nature 
of estimated parameters set. 

5 Diffusion model 

We start our analysis with the classical Black Scholes model, hereafter denoted BS, see Black 
and Scholes [3|. For the sake of simplicity, we take the money market account as the numeraire. 
Let ffii, J 71 , Pi) be the historical probability space and B the associated Brownian motion. 
Fix n £ M., o"o > and the interval [0, T\. The dynamics of the risky asset under the historical 
probability measure Pi is given by the following diffusion in accord with the model of Black 
and Scholes: 

(5) dSt = St fi dt + St o"o dB t . 

In this framework, the price of a European vanilla option is now standard (see for instance |24j). 
This model presents many advantages, in particular the pricing depends only on the volatility 
parameter ctq and we find closed forms for premiums and greeks of vanilla options. However, 
the BS model cannot reproduce the market price of call options for all strikes with the same 
volatility: this is called the smile effect. In order to take this phenomenon into account, we 
discuss two main extensions, the local and the stochastic volatility models, hereafter denoted 
LV, see for instance Dupire |13) . and SV, see Hull and White |20| . Heston |19| . Hagan et al 
|18] and Fouque et al |16j . In these two classes of models, the parameter o"o is replaced by a 
function a that depends on the time t, on the underlying S and, in the case of SV models, on a 
random source. For the sake of simplicity, we consider a local volatility model. The stochastic 
differential equation verified by the price of the underlying is 

(6) dS t = S t ndt + S t a{t,St)dB t . 

It is plain that this class of models is not the more general. But, your methodology can be 
easily applied to multidimensional diffusion, then the extension to SV models is possible under 
some hypotheses. Diffusions with jumps are excluded from the models treatable with your 
approach instead; for these, see for instance Cont and Tankov |10j . Henceforth, we denote by 
(Oi, J- , Pi) the probability space where the Brownian motion B is defined, by {J~t}t£[o,T\ the 
standard filtration generated by the Brownian motion B and by Ei the expectation under the 
probability Pi. 

Our goal is to analyze the sensitivities of the model given by the SDE ([6]) with respect to 
the uncertainty in the estimation of their parameters. Indeed, all models depend on certain 
parameters, generally a small number, that are related to the underlying, and this makes it 
possible to calibrate the model. When we select a given value for a parameter, using a calibration 
methodology, there is still some uncertainty on the true value of this parameter. This uncertainty 
can be estimated by using statistical methods, such as Fischer information, or by computing 
the sensitivity of the model calibration. We emphasize that both methods, statistical and 
calibrative, yield parameters with uncertainties and these uncertainties have a random character. 

Our goal is to analyze the impact of these uncertainties on the management and the hedging 
of a contingent claim from the point of view of a seller that tries to minimize his own risk. 

We propose to consider the model given by the SDE ([6]) with uncertainty on the parameter 
modeled by means of an error structure on the volatility function a(t, St). 

We make the following financial hypothesis: 

Assumption H-l (Asset evolution and uncertainty impact) 

1. the market follows the SDE (G|) with fixed but unknown parameters, i.e. the underlying 
follows the SDE |6j). The market is viable and complete, i.e. there are enough traded assets 
to guarantee that any contingent claim admits a hedging portfolio and there is a probability 
measure Q, equivalent to P, under which all discounted security prices are martingales; 



2. the option seller knows that the underlying follows an SDE such as §D\) but does not know 
the values of the parameters of the function a; 

3. the option seller has to estimate the parameters of his model: the uncertainty associated 
to such estimation is modeled by means of an error structure. As a consequence, the 
price and the greeks of the option are affected by uncertainty. We assume that the option 
seller makes his statistical estimations before to sell the option and he does not change his 
estimators during the time interval [0, T]. 

We assume that there exists a basis, denoted by {(fri(t, x)}j S n, of the space of volatility 
functions. The functions (f>i belong to C ' in (t, x), , {</>i}igN an d their derivatives are Lipschitz 
by a same constant and square integrable in x. For instance, the basis {4>i(t, x)}igN can be 
the spanning set of interpolating splines for the volatility function. We define the set L# of 
admissible volatility functions: 

Ll = I f{t, x) =^ j bi(/)i(t, x) such that (b\, b 2 , . . . ,b n , . . .) G f and f(t, x) > £ > V(t, x) > 

We equip the set with the distance a% induced by the sequence representation. 

We denote by cjj the coefficients of the series expansion of a, i.e. o~(t, x) = Y2% a i ^ii^i x )- We 
assume that the option seller has to estimate the coefficients a^, then the uncertainty is carried 
by the coefficients aj that are estimated by the estimators Ai, that is, the Ai are random 
variables. The estimated values of the coefficients a, are denoted by Ai, i.e. Ai are fixed values. 
We denote by <;(t, x) (resp. <f(£, x) ) the volatility estimator (resp. the estimated volatility), 
i.e. q(t, x) = ^ ii Ai(pi(t, x) (resp. ?(£, x) = X)i -^ <M*j x) ). We make now the following 
mathematical hypothesis: 

Assumption H-2 (Uncertainty on volatility) 

We assume o~(t,x) and (>(t,x) belong to L? a . For each random variable Ai, we define more- 
over an independent error structure (R, £>(R), Pi, B,j, r,j I . We denote by (Cl, T, P, B, TJ the 

product of all error structures (P, B(M), Pi, Bi, TA , that is again an error structure, see ^ 

chapter IV. We denote by {A, T>A) the related generator and its domain. We assume that the 
following hypotheses hold for each Ai. 

1. Ai G T>A, rL4j] and A[Ai] are known; 

2. the error structure (U, T, P, B, Tj admits a sharp operator denoted by (•) , in the sense 
that each error structure (R, B(R), Pi, Bi, Ti I admits a sharp operator denoted by (•)* 



and the sharp operator on the product error structure is the sum of the sharp operators of 
each sub-error structure. 

Finally, we assume the technical property 4-2. in Bouleau Jfi$ chapter V page 84- 

Property 4.2 in Bouleau [6] fix the probability measures Pi as a mixture of a measure \i 
absolutely continuous w.r.t. Lebesque one and a Dirac mass. Assuming the sum of the weights 
of the measure \x into the mixture is finite, then we can prove, see [6], that the probability 
product measure P = <g)Pj is absolutely continuous and under P only a finite number of the 
coefficients Ai in the representation are random variables. That means that only a finite, but 
random, number of terms contribute to the uncertainty. 

We briefly comment on the two previous assumptions. We have split them according to 
the nature of the assumption: the first one, IH-H is mainly financial and the second one is 
mainly mathematical. Assumption IH-21 sets the mathematical framework used in this paper; 
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point 4 and the final assumption are included in order to simplify proofs and to simplify some 
computations, but they can be relaxed. The hypotheses on the functions (f>i(t, x) are requisite 
to guarantee the existence and the uniqueness of the solution of SDE ©. The lower bound £ 
guarantees that the risk neutral probability Qi is well defined, in particular it is equivalent to 
Pi. 

The main financial hypotheses are contained in assumption IH-ll The first point states that 
the underlying follows an exact SDE without uncertainty. The uncertainty appears when the 
option seller attempts to determine the parameters of the SDE. This assumption is easy to 
understand in the financial framework and circumvents many problems in the mathematical 
framework. In particular, when the volatility is uncertain, the set of probability measures 
describing the whole class of possible probabilistic viewpoints is not dominated. In this case, 
the classical approach used in mathematical finance (see for instance Avellaneda [2]) cannot 
be followed; some important attempts to address this problem are Denis and Martini |llj and 
Denis and Kervarec |12| . 

When the option seller has set his parameters, i.e. his function <f(i, x), he has established 
his "risk neutral" probability Q>\, that is, the probability measure under which the diffusion 

(7) dxf ) = ft Xf ] dt + X ( p q (t, XP) dBt 



is a martingale, the related Brownian motion is denoted by W t . The probability Q\ exists and 
is unique since the risk premium is given by ft. We recall that the estimated volatility is a fixed 
function depending on time and on the underlying S, then it is known at any time given the 
value of the spot. Moreover, the estimated volatility <f is bounded from below by the positive 
constant £, then the Radon-Nikodym density is a martingale thanks to the Novikov criterion 
and the inequality 



Ei 



1 yV 

2 Jo ? 



exp {- J ^dt 



<Ei 



exp < — Todt 



2 Jo e 



< oo , 



Thanks to this probability, he can calculate the option price and the related hedging strategy. As 
a matter of fact, the replacement of the volatility function estimator q(t, x) with the estimated 
value of the volatility <f(i, x) causes an information loss. The natural question that arises is: 
What is the impact on the prices of the uncertainty on the volatility function <j(t, x) and it is 
possible to bound the related risk? 

We will analyze the problem in the following section 15.11 and we will propose a pricing 
principle in section [6) We shall also make the following assumption on the class of contingent 
claims analyzed in this paper. 

Assumption H-3 (Contingent claims) 

Let <I> be the payoff of a considered contingent claim. We assume that <I> depends only on St 
and belongs to C 2 . We also suppose that the two first derivatives of & are bounded. 

This hypothesis is needed in order to apply error theory techniques. However, we will 
show at the end of section [6] that this hypothesis can be overcome. Another point where the 
hypothesis about the boundedness of the first derivative of <3? is crucial is to prove that the gain 
process of the option seller is a Q^-martingale. 

5.1 Profit and loss process 

We consider that the option seller uses his data and proprietary information to determine the 
model parameters using an optimization procedure, see for instance Dupire |13| and Cont and 
Tankov |10] , By Assumption IH-lj a risk neutral-measure Qi exists and is unique. With this 
probability measure, the seller could define the fair price of the option and determine the hedging 
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strategy if he knew the true values of the parameters. Since he does not know these true values, 
he has to estimate them before to sell the option in order to propose a price for the option and 
to define an hedging strategy. 

The natural question is: how can he define this price taking into account the uncertainty 
into his estimated parameters? To do that, we recall that the wealth minus the liability of the 
option seller, that sells the option and hedges it, is the crucial stochastic process to define the 
price of the option, since the fair price and the right hedging strategy are such as the wealth 
minus the liability of the option seller is worth zero almost surely in a complete market without 
uncertainty. Therefore we study this wealth reduced by the related liability also know as the 
profit and loss process. 

We neglect, in this first work, the fact that the option holder can sell the option before 
maturity, so we assume that he holds the contingent claim until maturity. We also assume that 
all prices are denominated using the risk-free asset as numeraire. The profit and loss process at 
maturity of the option seller is given by the cost of the hedging strategy of the option, i.e. the 
option price in a financial model without uncertainty (see for instance Black and Scholes [1]), 
plus the value of the hedging portfolio minus the final payoff that the seller has to pay to the 
holder, so we have 



(8) 



PkL{T) = C(q, x, 0) + f A(?, St, t) dS t - $(S T ) 

Jo 



where C and A denote respectively the cost of the hedging strategy followed by the option 
seller, and its first derivative with respect to the underlying. Thanks to the asset pricing theory 
and in absence of parameter uncertainty, we have C(a, x, 0) = Ef 1 [$>(St)]- However, since the 
option seller does not known the true values of the parameters, the cost of hedging strategy is 
given by 



C(q, x, 0) = E' 



x®=x 



$ X 



-(?) 



where E K^ denotes the expectation under the probability space (Qi,J- , {•7 r t 1 }te[0,T]> ( 



X 



when the process X® starts at time with value x and follows the SDE (J7|), i.e. 



(9) 



dX. 



(C) 



xPt(t,x 



^) 



dW t 



(0 



dC 



dA 



under probability Q? . To simplify our notations, we denote by — — ——, - and r 

da{<t>i) da(4>i) 

derivatives of C and A with respect to the volatility in the direction of (j>, that is: 



the Gateaux 



< 10 > 8^-f' '• °» 



dA 
(11) tttttCS S t , t) 



da(<pi 



lim 






$ X 



-(?"+£ </>i) 



E" 



x, 



(?)_ 



$ X 



-(5) 






$ X 









$ x 






lim ■ 



dx 



dx 



x=S t 



The next lemma shows that the two previous Gateaux derivatives exist. 



Lemma 5.1 (Existence of the Gateaux derivatives) 



Under assumptions \H-^ and \H-^. the Gateaux derivatives defined by relations M0\) and Ml\) 
exist. Moreover, the higher order Gateaux derivatives — — , , N - — 7—^ and — — , , N - — 7—^ exist 



too. 



da{(t)i)da{^j) da(</)i)da(4>j) 
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Proof: We remark that A t is a Q^ martingale. Then, the process A t under Q^ has the 
same law as X t under Qi, where X t follows the SDE: 



dX 



(?) _ vK?) 



(?) 



xr? it, x?'\ dw i 



where W is a Qi Brownian motion. Therefore, we can write the difference quotient of C(?, x, 0) 
with respect to a variation of the volatility <f in the direction <j>i in equation (|10p in the following 
way 



sc _ 




'$ (xf** )' 


*0 =* 


[•( 


X?)] 


e 


Using: Assumption IH-31 we expand the function $ and it exists 


4 ? ' 0i,e) e" 


mm Ay , Xrp J, max 1 Ay , X T J 


5 


such that 













(12) 






E' ; - 



■(?, x, 0, e) 



^^ v(?.«^.e) + l<jj" (j{s,4>i,z) 



&IW)Y T 



Yn 



(?,4>i,t) 



T 



where the expectation is taken assuming the starting condition A - w 



(?+e «=XP=xand 



1- 



(?,</>i,<0 



W?+e</>i) W?) 



A' 



A 



We study the SDE satisfied by Y± 



(?>0i,<O. 



dY t 



(?.0i>e) 



A 



(?+e0*) 



[i,Af +e ^)+^f*^f +e ^ 



dWt 



-Af eft, xf> ) 



with K 



(?">&>«) 



ii 



0. The previous SDE can be written as 



dY 



(S,<f>i,e) 



eX 



+X 



(?+£</>;) 



t,xf +eM ) dW t 



t 



t[t,x 



(t+efc) 



i ^t 



t[t,x. 



(?) 



> Vi t 



dWt 



+Y f^k(t,xf ) )dw t . 



Assumption IH-21 guarantees enough regularity on <f to expand the second term. In particular, 
it exists 



J 



■(fj&.e) 



<E 



-(S+e<t>i) W?) 



min(Ar ^ ,Ar , max ^ ,^ 



K?+ £ 0O v<f) 



such that 



dF t 



(?>0i,e) 



-($+ £ <t>i) 



eXr^'^(t,X t 



(?+£&) 



rfW^ 



+ r 



(?,^,e) 



^Af)+Af +£ ^f (t,7f^ 



d^. 



Moreover, it is possible to choose J ' * " to be path-continuous thanks to the regularities of <f 
and the path-continuity of A. Then J is predictable. 
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We define the following processes: 

rt 
'o 

/ w\ r?x,j,.i r)c~ / (■? a.. c \\ 

dW K . 



H t Ixf+^frfsiXf^dWs 



Z t 

/o 



: r( S) x?)+xf +e «g( s ,7p^ 



We remark that -ff is an adapted process and Z is a continuous martingale, then we can apply 
the variation of constants method (see Protter |26] Theorem 53 section V.9) to prove that 

Y&Kd = eHt + e£ ( Z)t / £( Z )-1 (ifsrfZs _ ^ sd[Z) Z ] g ) ; 

Jo 

where £(Z) denotes the Doleans-Dade exponential of Z. In particular, we remark that Y^' oc 
e, then we can replace Y^ e into equation (|12p . Thanks to Assumption IH-31 the norm of 
the second term in equation (|12p is controlled by ce 2 , where c is a constant. Finally, thanks to 
the continuity of all processes, the difference quotient converge when e goes to zero, then the 
existence of the limit (fTU|) is proved. 

To prove the existence of the Gateaux derivative of A(q, St, t) defined by the limit (|lip . 
we have to study the SDE verified by the first derivative of &(X t ) with respect to the initial 
condition x. This SDE, the existence and the uniqueness of its solution and the closed forms of 
the solution can be proved adapting the proof of Theorem 39 section V.7 in |26j and is omitted. 
Then, the proof of the existence of the limit (jlip and of the second order Gateaux derivatives 
is similar to the previous one. 

□ 
We remark moreover that the integrability of A w.r.t. the asset S and, therefore, the fact 

r 

that the gain process / A(q,St,t)dSt is a Qi -martingale is guaranteed by the existence and 
Jo 

the uniqueness of the SDE solution verified by the first derivative of $(X t ) with respect to the 

initial condition x combined with the boundedness of the same first derivative, see hypothesis 

E31 

The previous result has a direct consequence: 

Lemma 5.2 (Gateaux derivatives as continuous maps) 

m „ dC dA d 2 C d 2 A 

The Gateaux derivatives ——. — r, — — ; — r, ——. — and as junction oj 

da(4>i) da(4>i) da{4>i)da{(t>j) da{cpi)da{(l)j) 

estimated volatility are continuous maps from the metric space (L^d 2 ) into the space of squared 

integrable function at the point <f. 

To prove this lemma, we remark that L 2 is an open set of the space of functions spanning by 
{</>i}igN equipped with the distance induced by the sequence representation. The target space 
is a Hausdorff space, then we have only to prove that the limits of our Gateaux derivatives 
depending on the function b(t, x) when b(t, x) approaches q are the same Gateaux derivatives 
evaluated on q. This result can be proved adapting the proof of lemma 15.11 and is omitted. 

The key remark about P&L equation ([8]) is that the price C and the portfolio A depend 
on the volatility process <f estimated by the trader. This estimated volatility is generally is 
different from market volatility a. Moreover, the estimated volatility is just the result of an 
estimation rule applied to the data set. The profit and lost process suffers from the uncertainty 
related to this estimation, since the estimated volatility q must be replaced by the estimator q. 
We concentrate our attention on the law of the P&L at maturity. In the absence of parameter 
uncertainty, the random variable P&L(T) is equal to zero almost surely and the option can be 
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hedged exactly. However, the exact hedging strategy cannot be implemented since the option 
seller does not know the parameters values aj but only the estimated A4. Therefore, P&iL{T) 
is a random variable and we have the following remark. 

Remark 5.1 (Random sources) 

The value of the profit and loss process at maturity is a random variable depending on two 
random sources: 

1. the stochastic "true" model (Oi, J- 1 , Pi) since the trader cannot use the correct hedging 
portfolio. 

2. the space (fi, J-, P), i.e. the stochastic process s, that depends on a random component 
unrelated to the Brownian motion Bf. 

As a consequence P&cL(T) is a random variable in the product space (f2, J-, P). 

Remark 5.2 (Role of the historical probability) 

The profit and loss process must be studied under the historical probability P = Pi x P. As a 
matter of fact, the risk neutral probability Q = Qi x P can be used if and only if all contingent 
claims are attainable. In our case, the option seller does not know the actual diffusion coefficients 
of the underlying, so the law of P&zL is not degenerate. The main impact is that the drift of 
the SDE {§[) plays a role and that the second term in the P&iL process J2J) is not a martingale. 
This fact complicates the computations in this paper. The important role of the drift fj, in asset 
pricing when the market is incomplete has already been emphasized in the literature, we mention 
for instance Karatzas et al. ]22j and Lyons f25\j . 

The two previous remarks show that the computation of the law of the PSzL process is not 
immediate. Moreover, the essential ingredient needed to determine the law of the P&L is the 
law of <j, and the law of this process is difficult to determine since it depends on the calibration 
methodology. We assumed in Assumption IH-2| that the option seller can estimate the mean 
value of the parameters and their variance, i.e. the two first moments of the law of q. 

Given the knowledge of the two first moments of the law of ?, we can at best estimate the 
two first moments of the law of P&L(T). This remark justifies our use of the error theory using 
Dirichlet forms, see section [3j As a consequence, we can define accurately the law of PSzL(T) 
with respect to the random source on the probability space Pi but we can just estimate the two 
first orders of the dependency with respect to P. In order to analyze the law of the PSzL process, 
it is sufficient to study the Pi-expectation on a class of regular test functions h(P&cL(T)) and 
the uncertainty on them using error theory. In practice we will compute the bias and the 
variance of 'Ki[h(PSzL(T))], where Ei denotes the expectation under the historical probability 
Pi. As a consequence of the dependence on s, Ei [h(PSzL(T))] is a random variable on the 
probability space (Q, J-, P). 

For sake of simplicity, we cut the dependency of P&L on T throughout the rest of the paper. 

Theorem 5.3 (Approximate law of the profit and loss process) 

Under assumptions \H-1\ \H-2\ and 177-ffl and for all test functions h belonging to C 2 with 
bounded derivatives, we have the following bias and variance: 

(13) A\Ei[h(P&L)]](a,?) = ^A 4 (1 V,?M[M^) 

i 

+^^(a,^T[A i ,A j }(A i ,A j ) 

(14) r[Ei[/i(P&L)]](,7,e) = X>i(<7, ?)*>, ^T^A^AuAj), 
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where 






Ei 
+Ei 



BC r T 8A 



h'(PkL) 



d 2 c 



da{^)i)da{(j)j) 



(?, x, 0) + 



d 2 A 



da(4>i)da( 



-(?, 5 t , i) dS< 



^"(P&L) <j — -(?, x, 0) + / 7r - TT (?, 5 S , s) dS s [> : 



5cj(0j 



<9cr(>;) 



*i(ff, ?) 



Ei 



x {^y(^'°) + r^) (? ' 5 - s) ^ 



Proof: For sake of simplicity, we neglect the dependency in T of P&L. The proof is split 
into two parts, we first prove the relation (|14|) and then (|13|) . Thanks to Assumption IH-21 the 
operator T admits a sharp operator. In particular, it exists a sharp operator denoted by (-W 
on each sub-error structure associated to each random variable A{ and the sharp operator (•)* 
on the global error structure is the sum of the sharp operators in each sub-error structure. 

We study the sharp of Ei [h (P&L)]. Thanks to the linearity of the sharp operator and the 
smoothness of the test function h we have 



(Ei[/i(P&L)]) # = Ei 



(15) 



E 



ti {P&L) (PkL)* 

ti (PbL) U7(s, x, 0) + J A(?, S t , t)dS t - $(St)J 



#' 



Using the linearity of the sharp, the term into braces can be rewritten as 

rT \ # 



(C(s,x,0)f+(J A(s, S t , t)dS t ) -mS T )) 



# 



are we remark that ($>(St)) = since the payoff is independent on the volatility estimated 
by the option seller. The integral is a linear operator defined by a L -limit, we proceed by 
approximation replacing the integral by a sum. We recall that the "true" diffusion S does not 
suffer from the uncertainty on estimated volatility. Using the linearity and the closedness of the 
sharp operator, we can then write 



J A(s,St,t)dS t J =J (A(?, St, t)f dS t , 



where we have used the first point of Assumption IH-ll see proposition V.8 page 83 in [6] for 
a more detailed analysis. Thanks to the expansion of the volatility (see Assumption IH-2] ) . the 
linearity of the sharp operator (see Proposition 13. ip . we have 



(16) 



C(C x, Of 

# 



£^(^0Mf 



I A(c;,St,t)dS 
Jo 



E- 4 



# 



T dA 



<9cr(c^ 



(?, St, i) dS t 



We multiply by h'(P!kL) and we take the Pi-expectation of the two right sides. Thanks to 
Assumption IH-21 Af is defined on a probability space (Q, x (7) distinct from (fii, F l , Pi), then 
we have 



E 



E 



ti{P&L)C{q, x, 0) # 
#■ 



ti(PkL) ( J A(q, S t , t) dS t 



i 



r)C. 



.4 



# 



T BA 
h'(PkL) I -—-( q ,St,t)dSt 

Off(^) 



.4 



# 
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then 



(E 1 [^(P&L)]) # =^^(cr, ? )Af. 



We conclude the proof on the relation for the quadratic error r [ E [h (P&L)] ] using the first 
property of the sharp operator (see Proposition [3TTJ, i.e. 






E 



44 



i.j 



and we find finally the relation (|14p evaluating the last random variable on the particular 
realization corresponding to observed data. 

The study of the bias is more complex, we first apply the closedness of the bias operator 
and the chain rule ([2]), then we find 



(17) A [Ei [h(PkL)}} = Ei [A[h(PkL)]\ = Ei 



h'(PkL) A [PkL] + -ti'(PkL) T [P&L] 



We study the two terms separately. We first consider the term depending on the variance. 
Thanks to the previous analysis, in particular relations (|16p . and the properties of the sharp 
operator, we have 

( r)C f T BA ~] 2 



'j 



da(4>i 



da(4>i) 



We multiply by h"(P&zL), we take the expectation under Pi and we find the third term in (1131) . 
We now study A[P&zL], we apply the linearity of the bias operator and we find 



A[PkL] =A[C(s, x, 0)]+A 



f A(s,S t ,t)dS t 
Jo 



AMS T )]. 



The last term is worth zero, since the final payoff is completely defined by St and does not 
depend on the volatility estimated by the option seller, see Assumptions IH-ll and IH-31 Thanks 
to the same assumptions and the closedness of the bias operator A, we have 



A 



[ A(<;,St,t)dSt 

Jo 



T 



A [A( ? , S u t)] dS t . 



Thanks to the same argument used for the study of the variance, we can take the Gateaux- 
derivatives of C and A with respect to a variation of the volatility along the component fa and 
using the chain rule ([2]), we have 



A[C(<;,x,0)] = Yl 



dC 



4£ 



fax, 0)A[Ai 

d 2 c 



I 



, A^x,0)T[A i ,A j ] 

[' fT flA 

A[A(,,S t ,t)]dS t = E y d^^> St > t)dSt 



A[A] 



4E 



i -j 



T 



d 2 A 



da(4>i)da{<j)j 



■(s,S t , t) dS t 



VlAuAA. 



We multiply by h'(PkL), we take the Pi-expectation and we find the two first terms in (fT3 
evaluating this random variable on the particular realization corresponding to observed data, 
thanks to the continuity at the point (A±, . . . Aj-, ■ ■ •)■ 
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□ 

We remark that the approximate law of the profit and loss process depends both on the 
true volatility a(t, x) and on the estimated one q(t, x). Two natural problems arises. First of 
all, the true volatility a is unknown. Moreover, the volatility estimator q is a random variable, 
hence the bias and the variance of Mi[h(PSzL)], given by Theorem 15. 3| are conditional moments 
known given the expected volatility ?. It is important to remark that ¥.i[h(PkL)] remains a 
random variable on the probability space f2. We explicit this fact indicating the dependency on 
? throughout the rest of the paper. 

In a particular case, we can partially overcome the problem of the dependence on a. We 
have the following corollary 

Corollary 5.4 (Bias and Variance of the expected value of PkV) 

Under assumptions \H-l\ \H-2 and \H-3i we have the following bias and variance: 



(18) 



(19) 
where 



i 



>;l 



r pEi [PkL]] (a, ?) = J^ **( CT > *)%(*' 3 r ^< A i](Ai,K 
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A?V,3 



A&V.9 



M<r, ?) 



dA 



dC f T 

— — — -(?, x,Q)+fi E ■ 
da((j)i) J [da ((pi) 

d 2 C f T 

da((j)i)da(4>j) J 

dC f T 

— -t(^ x, 0)+n / Ei — — 
oa(4>i) J [oa((pi) 



(?, S t , t)S t 



dt 



d 2 A 



d[°(&)] 2 
dA 



(?, S t , t)S t 



dt 



(q, S t , t)S t 



dt 



Proof: The proof is a direct consequence of the Theorem 15.31 where we have calculate the 
expectation using the SDE ©, we omit the details. 



□ 

In particular we remark that A t (a, <f), A i (a, <f) and ^i(a, <f) depend on a by means of 
the portfolio term. Moreover, assuming /i = 0, the dependence on a vanishes. In other words, 
the bias and the variance of the expected value of the option seller wealth are independent on 
the true volatility in a market without risk premium. 

In the general case, we can approximate the equations (|13p and (fT^|) using the known value 
<f instead of the unknown a. Under this approximation, we have the following corollary. 



Corollary 5.5 (Approximated tails of the law of the P&.L) 

Under the hypotheses of Theorem 15.51 and assuming the approximation a = q, the random 
variable Ei [h(PSzL)] (q), defined on probability space (£1, J 7 , ¥), has the following approximated 
tails: 



(20) P Ei [h(PkL)] (?) - /i(0) - A [Ei [h(P&L)\] (?) > k yjT [Ei [h (PkL)]} (? 



< 



1 

1 + k 2 
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where 
(21) 



i 

+l 1 £ t A?](t,t)T[A i ,A j ](A\,A i 



'■j 



(22) 
with 



T [Ex [h (PkL)}} (?) = ]T ^(c, ?)^-(?, ?) T[Ai, Aj](Ai, Aj), 



ij 



A - 1,(f ' f) - ft ' (0, fe ( " >0)+ ''i El 



xp^ i? ,xr,t) 



^) 



da{(pi) 



,( 2 ) 



A$(S?) = fc'(0) 



9 2 C 



da{4> i )da(4) j 



-(?,x, 0) + n [ Ei 



9 2 A 



r// 



' cV(<^)<M<A 



i?,xr,t) 



(?) 



■j> 



r// 



x W ( " o) W * 



X® dA (c X® s) 



ds 

ds 



+h"(0) / E x 
Jo 

8(1 

*,(?,?) = ti(0) 



{xP} 



2 ^ (? X® S ) ^ (r X ^ s) 



da(cf>i) 



(?, x,0)+n Ei 

JO 



X (f) 



dA 



da((t)j] 
- X® 



ds 



ds . 



TTie right-hand member in relation A20\) can be replaced by the complementary error function 
assuming the Gaussian approximation. 

Proof: We first make an useful observation: In a complete market without uncertainty, the 
price C and the portfolio A are defined in such way that the P&zL process is worth zero almost 
surely. The changing of a for the estimated value <f into the definitions of A\ and ^u\ entail 
that the PSzL is substituted by the approximated process 



-(?) 



^ 



P&L = C(?, x,0)+ A(q, Xf>, t) dXl - ^(X^>). 
Jo 



Hence, by definition of C and A, we have that P&iL = Pi-a.s.. We know the bias and the 
variance of Ei [h(PSzL)] thanks to Theorem l5.3l Then, we approximate them using a = <f and we 
find the approximated bias .4[E[/i(P&:L)]](?) and the approximated variance T [E [h (P&zL)]] (?) 
thanks to an easy calculation using the SDE (j9]). Finally, we apply Chebyshev's inequality, see 
Corollary 13.31 and we find the relation (|20p thanks to the previous remark. 

□ 

The last corollary gives us the approximated tails of the law of the P&L, in particular it 
permits to estimated the probability of a rare event thanks to formula (|20p . In the next section, 
we will use these results in order to define the price of the contingent claim. 



6 Option Pricing 

In order to interpret this result in financial terms, we consider that the option seller is aware 
of the presence of errors in his procedure for estimating the volatility a and wants to take this 
into account. Since the option seller does not control the errors, the risk related to the space 
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(fi, J 7 , P), i.e. to the uncertainty on the parameters, cannot be hedged. Indeed, if we compute 
the super-hedging price of a contingent claim, we find a very high buy-price and a very small 
sell-price, i.e. the bid-ask spread becomes too large compared with the market spreads. We 
also remark that the probability that the random variable E[(P&L)] takes values far from the 
mean is very small and becomes negligible if this distance is big compared with the standard 
deviation -y/r[E[(P&L)]]. Therefore, we introduce the following principle to define the price of 
a contingent claim under this uncertainty model. 

Principle P-l (Asset pricing principle under uncertainty) 

The option seller fixes a risk tolerance a < 0.5 and accepts to sell the option at any price 
Fsell(s)> J^o-measurable, such that 

(23) P { Ex [F, rf ,(?) - C(s) + P&L(?)] < 0} < a 

where C(?) denotes the cost of the hedging strategy, that is the theoretical price of the option 
without uncertainty on the volatility, i.e. C(?) = E^i[$(X^ )]. 

Before applying this principle to our analysis, we discuss its financial implications. First of 
all, we remark that -F se ii(<f) and C(?) are J-q measurable, then the relation ([23]) can be written 
as 

(24) P {FseiiOf) - C® + Ei [PkHs)\ < 0} < a . 

Moreover, the price .F se n(<f) depends on the estimated volatility <f. The second remark is that 
this principle does not give a single price but a half-line of possible selling prices. Indeed if X 
is a selling price and Y > X then 

P {Y - C(?) + Ei [PbL(q)] < 0} < P {X - C<£) + Ei [PkL(q)] < 0} < a, 

so Y is also a selling price. We call ask price, denoted i ? a sk(?)> the infimum of the set of all 
acceptable selling prices. We also remark that the classical option price in the model without 
uncertainty is a sell price and, in particular, is the ask price, since the wealth of an option seller 
is worth zero and the probability space P is trivial, in a complete market without uncertainty. 

We discuss the event {Ei [X — C(?) + PSzL(q)] < 0}. On this event, if the option seller has 
sold the option at the price X, then he will lose money in mean at maturity, since he has to 
pay C(?) to buy the hedging portfolio and the noise on this hedging strategy conducts him to 
lose money in mean. However, the probability of this event is smaller than a. 

Therefore, the asset pricing principle says that the option seller accepts to sell the option 
at a price X if this price is high enough to guarantee that he loses money with a probability 
smaller than a. The choice of the parameter a depends on the risk aversion of the option seller. 
However, if he is too risk-adverse, then he will advertise the contingent claim at too high a 
price, and potential buyers will find other traders that offer the same option at a lower price. 

We remark that there is a likeness between this principle and the VaR of the trader portfolio. 
Indeed, a can be interpreted as the probability to "lose money on the contract", which is the 
definition of the Value-at-Risk. We also remark that the final wealth of the trader is negative 
in mean on the event {Ei [X — C(?) + P&iL(q)\ < 0}; this event is the most onerous then we 
name the related risk the "risk of the first kind". The "risk of the second kind" for the trader 
is to propose a price too high such that potential buyers do not buy the option, i.e. "lose the 
opportunity to make money" with the contract. In a future paper, we will study the problem 
of the optimal proposed price as an equilibrium between the two previous kinds of risk. We 
assume here that the risk tolerance a is fixed and that the option seller chooses the smallest 
price consistent with his risk tolerance. 

We also remark that the Principle IP- 1 1 also defines the purchase price. As a matter of fact, if 
a trader accepts to buy an option, he has to take a negative position on the hedging portfolio to 
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cover it, i.e. he has to follow the opposite of the hedging strategy. In accord with the Principle 
IP- 11 the option buyer accepts to buy the option at any price i ? buy(?) such that 

P{-F buy (e) + C(?) - Ei [P&L(0] < 0} < a, 

that is the relation (|23p . where we have changed all signs. The relation above can be rewritten 
as 

SUiWOf) - C(s) + E x [PbL(s)] > 0} < a. 

The previous relation and the relation (|23|) lead us to the following remark: 



Remark 6.1 (Bid- Ask Spread) 

If all option traders on the market are risk adverse, then the best purchase price, denoted by 
Fbid(s)> an d the best seller price -F as /c(<f) are distinct. 

This remark results from the fact that if the traders are risk adverse, then the risk tolerance 
must be smaller that 0.5, so F b id(?) < C(?) + E [Ei[P<kL(V)]] < Fask(<r)- 

The previous principle and Theorem 15.31 have the following immediate consequence: 

Proposition 6.1 (Option prices) 

We assume the hypotheses of Theorem 15.31 the approximation a = q~ and that the option 
seller follows the principle \P-ll Then he accepts to sell the option at any price higher than 



Fask® = C(?) + A pEi[P&L]] (?) + y/r\Ei[PbL]](3M 



1-Q 



where N\- a is the (1 — a)-quantile of the reduced normal law, or the function . h^— given by 

Chebychev's inequality ^ in the conservative case. Likewise, the option buyer accepts to buy 
the option at any price lower than 



Fbid® = C(s) + A [Ei[J>&L]] (?) + vTpEi[P&L]](?)M, 

We remark some symmetry in the two previous prices, since M a +N\- a = 0. Therefore, we 
have 

(25) F mid (s)= E * ski ' )+ 2 Fhidi < ) = C® + A\Ki[PkLm 

(26) Spread(f) = F ask (?) - F bid (?) = 2 y /r\E 1 [PkL]]{?)M 1 - a 

where i ? m id(?) denotes the mid price, i.e. the average of the bid and ask prices, and Spread(<f) 
denotes the bid-ask spread, i.e. the difference between the ask and the bid prices. We emphasize 
that with our model, we can reproduce a bid-ask spread and we can associate its width to the 
trader's risk aversion, i.e. the probability a, and the volatility uncertainty, i.e. the term 
y/T [EifP&L]] (<f). Another interesting point it that the mid-price does not depend on the risk 
tolerance a. 

We conclude this section with a remark on the Assumption IH-31 

Remark 6.2 (General Payoff) 

It is clear that the Assumption \H-3\ is quite restrictive, since call and put options do not 
verify it, for instance. However, we use the hypothesis that the payoff is regular only to assure 
the existence of the Gateaux derivatives of the price and the hedging portfolio, see equations 
H0\) . \11\) and Lemma \5.1\ if we can prove that the Gateaux derivatives exist, then Theorem \5.3\ 
remains valid and we can define the price for an option according to the Principle \P-1\ even if 
the payoff is not regular. 
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7 Example: log-normal diffusion 

In this section, we give an example of the previous results. In particular, we consider the log- 
normal diffusion, i.e. the Black Scholes model. The underlying follows the SDE ([5]). In this case, 
the volatility is a real number, so we replace it by the same number multiplied by the constant 
function. We concentrate our analysis on call options. It is clear that a call option does not 
verify the Assumption IH-31 but we have a closed form for all greeks in the case of a log- normal 
diffusion and we can check that the Black Scholes pricing formula is C 2 at each time t strictly 
smaller than the maturity T and the vega, i.e. the derivative w.r.t. the volatility, vanishes 
when t goes to T. These properties guarantee that Theorem 15.31 remains valid even without 
Assumption IH-31 see Remark 16.21 We assume that the drift \x vanishes under the historical 
probability measure Pi, in order to simplify our numerical computation and to overcome the 
problem of the dependence on ctq , see section I5.1) in particular Corollary 15.41 We denote by ft 
(resp. ft) the volatility estimator (resp. the volatility estimated by the option seller). For that 
we consider an error structure of type Ornstein-Uhlenbeck, see Definition 13.21 
In this case, the Theorem 15.31 has the following corollary. 

Corollary 7.1 (Bid and Ask prices with log- normal diffusion) 

If the underlying follows the Black Scholes SDE f5j) without drift, then the bias and the 
variance of the profit and loss process that hedges a call option verify the following equations: 

(27) A[C(x,K,T,, W ) = x ^Vt\a[,o}(9>) + ^ r [«,](*)) 

V2vr [ 2 a J 

-d 2 

(28) T[C(x, K, T, ft)](ft) = x 2 ^-TT[ % }(%) 

where 

lnx-lnisT+fr _ r - 

d\ = -j= — - — , d 2 = di- ft VT, 

ft VT 

and C(x, K, T, ft) denotes the price of the call with strike K, maturity T when the underlying 
has value x and the trader known the volatility estimator ft. Moreover we have the following 
bid and ask prices: 

-id? 

(29) C ask {x, K, T, ft) = xQ{d x ) - KQ(d 2 ) +x^fT^= jf [<&] (%Wi-a 

V 2n 



+x VT e -^ {a [ ft ] (%) + ^ r [ft] (ft)) 

V2vr [ 2o- J 



-id 2 
e 2 a i 



(30) c m ( x , k, t, ft) = x g{d x ) - k g{d 2 ) + x Vt —=r y/r [«,] (? )at q 

V Z7T 



+x VT e -j= [a [«] (a) + ^ r [ft] (ft)} 



where Q denotes the cumulative distribution function of the reduced Gaussian law and J\f a is 
defined in Proposition \6.1[ 



Proof: We know that the cost of hedging strategy -F(ft) and the portfolio A(ft) of a call 
on Black and Scholes model are given by 

C(x,K,TSo) = C(^) = xg(d l )-KG(d 2 ) 

dC 

A(x, K, T, ft) = _(ft) = g(d 1 ). 
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We can easily calculate the two first derivatives of F(%, x) and A (5), x) with respect to 
the volatility %, then the two equations (|29p and (|30p come directly from a computation of 
A[C(x, K, T, ?0 )](?o) and T[C(x, K, T, ?0 )](?o) using Corollary El 

D 

7.1 Analysis of the impact of parameter uncertainty 

We now analyze the correction on the pricing formula due to the presence of an uncertainty 
on the volatility parameter. We have the following corollaries that can be proved with simple 
computations. We analyze only the mid-price, i.e. the average price of the bid and ask prices. 
We also use the notation 

Corollary 7.2 (Delta and Gamma correction) 

The bias on the call price due to the uncertainty on volatility verifies 

rsn dA[c] - dl A[c] x e ~ ¥i ui +,nrr?i 

Moreover, the bias and its first derivative with respect to K are positive at the money ifr^ (5)) > 
j^q T and are negative when the stricktly inequality is reversed. Finally, the bias is convex in 
K at the money if and only if 

Corollary 7.3 (Time evolution of the correction at the money) 

The third cross derivative of the bias, twice with respect to the strike K and once with respect 

■ ■ d 3 A[C] . ...,,,., 

to the maturity, i.e. , is positive if and only if 

s l qrJgr-tf + 128 

r [%>> 4 16 + ^T 2 

The two previous corollaries have an interesting consequence. 

Proposition 7.4 (Smile on implied volatility) 

There exists an interval ]a, b[ such that for all r^ s (^Q) G]o, b[, the implied BS volatility 
recognized from the average price F m id given by Corollary 17.11 is a convex function around the 
money. Furthermore, the second derivative of the implied volatility with respect to the strike is 
a decreasing function of the maturity T. That is, the implied volatility exhibit the smile effect 
and this effect is stronger for short maturities than long one. 

Proof: We start fixing by rf (%) = j?qT. Thanks to Corollary 17.21 we have that 
A[C] = —fix = at the money. In addiction, the same corollary assures us that the bias 
A[C] is strictly convex around the money thanks the continuity of the derivatives. For the 
same value of r^ s (^o) we easily check that the cross derivative of the bias two times w.r.t. K 
and one time w.r.t. T is negative at the money. We argue that the call price at the money is 
equal to the Black-Scholes one whereas around the money the call price in our model is bigger 
than the Black-Scholes one. Giving that the Black-Scholes formula is used to find the implied 
volatility, we conclude that the implied volatility is strictly convex around the money and that 
this convexity is more marked for short maturities. We now remark that the implied volatility 
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belongs to C°° as function of the strike K since it is a composition of infinite differentiable 
function. We conclude that, for each sufficiently small neighborhood of the money, it exists an 
interval ]a, b[, including -^ T, such that the implied volatility remains strictly convex and with 
convexity decreasing with the maturity. 

□ 

8 Example: Constant Elasticity Volatility Model 

In this section, we present a second example. We consider the CEV model, see Cox and Ross 
[9]. The volatility function is equal to 0"5f~ where a is a positive constant and 6 is constant 
and belongs to (0, 1). As in the previous example, we suppose for sake of simplicity that the 
risky asset is a martingale under historical probability Pi. We assume then that St is bounded 
from below by a positive constant. The diffusion is then 

dS t = o S^dWt 

This diffusion admits a closed form solution via squared Bessel process theory, see for instance 
Jeanblanc et al. |21j . The previous SDE can be rewritten in the following form: 

dS t = e loga+l31ogSt dW t 

In this form, it is easy to adapt our previous result to CEV diffusion. In our case, we consider 
that the couple of parameters (a, B) is fixed but unknown. The option seller knows only the 
couple of estimators (q,B) and the estimated values (£", B). We fix an error structure for the 
couple of the parameters: 



(R+x(0,l), 



x(o,i)),p,d,r) 



The explicit form of the probability measure P and of the carre du champ operator are not 
compulsory since we need only to compute the bias and the variance at the point (q,B). We 
now show how to compute in our framework the option price in a particular case, that is a 
regular but non-linear payoff. 



8.1 Variance and bias of the diffusion process 

We consider the diffusion process from the point of view of the option seller: 

B 



(32) 



d*P> = q 



n 



(S,B) 



dWt 



We have the first two results about the sharp and the variance of X^' > : 

Proposition 8.1 (Sharp of the process X^' B ^) 
The sharp of X^' ' verifies the following SDE 



(33) d(x^)* = qB(x^ B ))*[x[ 



^,-B) 



B-l 



dW t + 



XI 



fc-B) 



t \ ) t 

Moreover, the previous solution admits the following closed form: 



Li 



? # + B #A dWt . 



(34) (X^ 



# 



? # + b#A M t 



X^ M~ l dW s -qB / X^ 
- 1 Jo L 



2B-1 



M7 l ds 



where M t = £ \ qB 




X k,B) 



B-l 



aw* 
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Proof: We omit the dependency of X^ ,B > on the couple (<j, B) for simplicity's sake. We 
start discretizing the SDE ([32]) using Euler's scheme. Let II n = {0 = t% < f{ < . . . < t n n = T} 
be a subdivision on the interval [0, T] for each integer n with the discretization step denoted by 
5 n = supj(tj + i — ti). We have 

f x? = x Vte [#,*?) 

Xfy = X t n n_ i +s(x r t 'n_) (w t n-W t n\ Vi = l,...n 
X} = xf l " 1 Vfc€[i?,#u). 

We now apply the sharp operator to our discretized process X n . On the first interval [£q,£"), 
the process X n does not depend on the couple (^, .B) then X^ = for any t € [£q , £"). At time 
ti, we have 

# / „ \# / / „ \ B ' ^* 



thanks to the linearity of the sharp. The second term of the previous relation is the sharp of a 
function depending both on the couple (?,-£>) and on the process X n . In order to compute it 
we apply the second property of sharp operator, see proposition 13.11 and proposition V.8 page 
83 in [6] for a more detailed analysis. We have for any regular function f 



(/ M.*i,)) # = §: (^.*i> # +;S: (**.*&-,) B *+H: M.*i.) K, 



where ^- denotes the partial derivatives w.r.t. the jth variable. In our case we find 

+,(x%_ i ) B (w t? -W^ 1 )B# 

+ SB (X% Y ^ (W t n - W t n ] (XVU 



L *?_J V vt i - "t! %, 



Then, we have 

x|) # = (x^y+^iy^-w^),* 

K fe_J B (W t « - Wt ?1 ) B* 

■KB (-X^J 5 " 1 (W t? - W^) (X|_^ # 

and (X t n ) # = (X$ J for any t E [£™,£™ +1 ). We remark that {X n )f is a Picard sequence that 



verifies usual conditions. As a consequence it exists an unique process (Xyf defined into the 
probability space P x P that is the limit of the sequence (X™y when 8 n goes to 0. This limit 
is the sharp of the process Xt thanks to the closedness of the sharp operator (•) . Moreover, 
{Xyf verifies SDE (j33|) Finally, (|33p is an affine SDE, then we can apply the method of constant 
variation, see for instance [26], to find a closed form solution. It is easy to check that equation 
flM} is the solution of SDE (jggjl. 

□ 



Corollary 8.2 (Variance of the process X^' B >) 
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The carre du champ operator acts on X^ ,B > as 



(35) 



X, 



(S,B) 



(S,B) 



Mf 



-* *■ e 



M^dWs - qB 






2B-1 



M7 1 ds 



Proof: Thank to the first property of sharp operator, see proposition 13, II we have that 



* 



(S.-B) 



E 



xw 



# 



Using the closed form (|34|) and recalling that only ? * and -B* depend on P, we have that 



r 



X 



(C,B) 



Mr 



xE 



UO 






B 



Af^dPVa - ?5 



f r 



2B-1 



X^ B ) M." 1 ^ 



The last expectation is easy since E (<> ) = r[<j] and E (-B*) = T[-B] by definition of 

sharp. Finally, thanks to the continuity of the carre du champ operator V at points if and B 
and using lemma [531 we find equation fj35[> . 



□ 



Corollary 8.3 (Covariance of the process X^' B ') 
We have the two following covariances 



(36) 



?,X 



fc-B) 



>^ v t 



(?,5) = M 4 



t r 



B 



X^ M-^dWs - ?B / xi'' B) M7 x ds 



t r 



2B-1 



(rW(?) + ?r[ ?( fl](?,5)} 



(37) 



B,X 



(f,B) 



(e,£) = M t 



L 



x(?' s ) 



Mr^Ws - ?B 



o L 



Xf' fi ) 



2B-1 



m: 1 ^ 



x{«T[B](B)+r[ f ,J3](?,B)} 
The proof is similar to the previous one. 

Proposition 8.4 (Bias of the process X^' B ^) 
The bias of X^> > verifies the following SDE 



dA 



X ^ B ) 



sB 



X, 



(C-B) 



B-l 



A 



x^ B ) 



dW t + 



X^ B) (A[<;] + <;A[B])dWt 



(38) 



+-,B(B - 1) 



X 



k,B) 



B-2 



+ 



X 



fc-B) 



B-l 



sr 



r 



X (^ s ) 






dW t + 



X 



(f.B) 



B 



r[c,B]dWi 



s,xfc fi ) 



dWi 



1 

+ 2 ? 



'(«>*) 



B 
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Moreover, the previous solution admits the following closed form: 



(39) A [X t 
where 

v(s,X s ,^,B 



(?■*) 



<S,B) = M t 



t - /"* r 

Z /(s,X S) c,S)M s - 1 dW s + / d zv(.,X,f, J B),M- 1 

JO L 



xf> s ) 



B 



1+^? 



+-tB{B - 1) 



+ 



xF' B) 



B-l 






sr 



B-2 

r 



x 



?,B) 



X (?,B 

?,£) + 



B 



r[B](s) 



-* *- c 



(c ;J B) + er 



BjM 



r[?,B](?-,B) 
M)} 



Proof: The proof is based on the same arguments used to prove proposition 18.11 We 
discretize the SDE (|32|) using the same Euler's scheme. Then, we apply the bias operator A at 
the discretized process X™. The chain rule ((2J) verified by the generator depends both on the 
bias and the variance. Using the chain rule and the linearity we have 



A 






A 






8 ( Wt n - W tU ) Aft + ? fey (W t n - W tJLl ) -4[B] 



+?B ( X£, . ) [W&-Wt».)A 



X fJ + ^ B ( X I- 1 ) B ( w % ~ W V-i) r ^ 



+^B(B-l)(x%_ i Y 



-KBlX, 



+B[X%, 



B-l 



B-l 



W t v- — Wf. 



W t v- — Wtv 



B , X+ri 
C, X"n 






+ x& 



B 



Wtn-w^jr^.c] 



We recall that T is a closed operator, then T [X t n ] (resp. T [B, X t n ], T [?, X™]) converges to T [X t ] 
(resp. r [B, Xt], T [<;, Xt]). Then, we have construct a Picard sequence A [X"] that verifies usual 
conditions. As a consequence, it exists an unique process A [X] t defined into the probability 
space P that is the limit of the sequence A [X™] when 6 n goes to 0. This limit is the bias of the 
process Xt thanks to the closedness of the bias operator A. Moreover, A [X] t verifies SDE (f38|) 
that is affine w.r.t. .ALX]^. Applying Ito formula it is easy to check that equation (|39p is the 
solution of SDE (|38p . where we have used the continuity of A and T at point (q, B) and lemma 
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Finally, we can resume some interesting results about the variance and the bias of the process 
X\ \ ' , that can be deduced easily from the previous results. 

Proposition 8.5 (Properties of the variance and the bias of X t ' ) 

We have the following properties: 

X t (<f, B) is a non-negative local-sub-martingale taking value at time 0. 



1. r 

2. r 

3. A 



a, 



teA) 



1 y ^t 



(<f, B) and V 



A,X t 



isA) 



($,B) are local-martingale taking value at time 0. 



X 



isA) 



(<f, B) is a local martingale. 
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8.2 Financial aspects 



The easiest example is an option that pays S T at maturity. The price without uncertainty can 
be easily computed using the exact law of St by means of squared Bessel process approach or 
applying Monte-Carlo techniques thanks to Euler discretization. In both cases, we have the 
(approximated) law of the process S for any time t G [0, T\. Suppose now that the option seller 
does known the true value of the couple (a, f3) but only the estimated value (<f, B) and the 

estimators (q,B). The option seller has to work with his estimated process X\ ' . Clearly the 
hedging cost can be estimated using the exact law given by Bessel process theory or thanks to 
Euler discretization. In both cases, it is possible to estimate the bias and the variance-covariance 
of X\ ' depending on the uncertainty on (<j, B). 



Exact law. Given the exact law of X+ ' , we can compute the joint law of T 



X 



fcs) 



-fcB) 



(?,B), r 



B,X t 



fcs) 



(?, B) and A 



X, 



fc-B) 



(?, B) using equations (f35]j . ([36 



3 



and (J39D- 

Euler scheme. We recall that the operators carre du champ V and generator A are closed. 



Then we can define an Euler scheme to find the approximated joint law of V 



X 



fc-B) 



&B) 



S,X, 



(*,b) 



(?,s),r 



B,X, 



(*,b) 



(<f, B). Moreover, we have used the same 



(?,5)and„4 X^' B) 

approach to find the SDE verified by the sharp and the bias. Then, we can use a unique 
partition II n of the interval [0, T] in order to estimate the hedging cost and the four processes 



X 



(s,B) 



(t,B),T 



S,X; 



(s,B) 



AX 



(?,S) 



(<T, B) and A 



We now analyze the bias and the variance of the square of X t 

Proposition 8.6 (Variance and bias of (X^ J ) 

We have the two following relations 



X 



(«,-B) 



foB) 



(S,B). 
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(S,B) 
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2X 



(S,B) 
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(?.*) 



(?.B) 



<S,B) 



(t,B) + T 
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(?.-B) 



(S,B) 



The proof is a direct application of the two chain rules ([I]) and ([2]) verified by operators T 
and A. 

We can now conclude with the price of the power option in the case of uncertainty. 

Proposition 8.7 (Bid and ask prices of the power option) 

We have the following bid and ask prices for the power option in accord with Principle IP- 11 



(42) 
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where N a are defined in Proposition \6.1\ 

The proposition is a direct consequence of Proposition I6.1 



M a 



(S,B) 
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